Under consideration for publication in J. Fluid Mech. 



1 



Reentrant Hexagons in non-Boussinesq 

Convection 

By Santiago Madrugaj, Hermann Rieckef, and Werner Pesch$ 

fDepartment of Engineering Sciences and Applied Mathematics, Northwestern University, 

Evanston, IL 60208, USA 
JPhysikalisches Institut, Universitat Bayreuth, D-95440 Bayreuth, Germany 

(Received 8th February 2008) 

While non-Boussinesq hexagonal convection patterns are well known to be stable close 
to threshold (i.e. for Rayleigh numbers R « R c ), it has often been assumed that they 
are always unstable to rolls already for slightly higher Rayleigh numbers. Using the in- 
compressible Navier-Stokes equations for parameters corresponding to water as a working 
fluid, we perform full numerical stability analyses of hexagons in the strongly nonlinear 
regime (e = R— R c /R c — 0(1)). We find 'reentrant' behavior of the hexagons, i.e. as e 
is increased they can lose and regain stability. This can occur for values of e as low as 
e = 0.2. We identify two factors contributing to the reentrance: i) the hexagons can make 
contact with a hexagon attract qr that has been identified recent l y in the nonlinear regim e 
even in Boussinesq convection l)Assenheimer fc Steinberg! 1 19961) ; IClever fc Bnssel ( 1996^1 ) 
and ii) the non-Boussinesq effects increase with e. Using direct simulations for circular 
containers we show that the reentrant hexagons can prevail even for side-wall conditions 
that favor convection in the form of the competing stable rolls. For sufficiently strong 
non-Boussinesq effects hexagons become stable even over the whole e-range considered, 
< e < 1.5. 
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1 Introduction 

Rayleigh-Benard convection has served as an excellent paradigm for studying systems 
that spontaneously form spatial or spatio-temporal patterns. In recent years exciting re- 
sults have been obtained for the stability and dynamics of structures that are connected 
with roll convection. To mention are, in particular, the spiral-defect chaos obtained in 
convection of fluids with low Prandtl number and domain ch aos driven by the Kuppers - 
Lortz instability in rotating systems. For a recent review see iBodenschatz et all (2000). 
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In most of these investigations great care has been taken to keep the experimental sys- 
tems close to the regime in which the Oberbeck-Boussinesq approximation is valid by 
minimizing the dependence of the fluid parameter on the temperature in order to avoid 
the appearance of cellular or hexagonal structures. 

It has long been recognized that variations of the fluid parameters with the temper- 
ature, i.e. non-Oberbeck-Boussinesq (NOB) effects, break an up-down symmetry and 
therefore introduce otherwise prohibited mode interactions. The best-studied case is the 
resonant triad interaction among the three fundamental Fourier modes whose wavevectors 
form a hexagonal pattern. It renders the primary bifurcation to the hexagons transcrit- 
ical with the con sequen c e that the hexagons are preferred over rolls in the immediate 
vicinity of onset l|Bussel l)l967|) ). According to leading-order weakly nonlinear analysis, 
hexagons typically become unstable to rolls further above threshold where the ampli- 
tudes are larger and the reson ant-tr i ad int e raction loses significan c e compared t o inter- 
action s involving four mo d es JPalml lll96f I: ISegel fc Stuartl l|l962|) : ISegel l|l965l) : iBussel 
1119671) : IPalm et all l|l967h : iDavis fc Segel (jl££$). This scenario of a transition from 
hexagons to rolls has been confirmed in quite a num b er of experimental investigations 
" Somerscales fc DougheTtvI l 197C k lDubois et al\ lll97^:lRichterfl)l978(l : rWalden fc Ahlerd 
198l|) : iBodenschatz et all 1 1991 ): Pampaloni et all l)l992f) ) and it has commonly been 



assumed that in convection hexagonal patterns are confined to the regime close to onset. 
Presumably due to this assumption, hexagonal convection in strongly nonlinear non- 
Boussinesq convection has not been investigated much. 

Recently, however, two experiments showed that even in the strongly nonlinear regime 
stable hexagonal convection patterns can be observed. In one experiment hexagons were 
observed at relatively high Rayleigh numbers (e = (R — R c )/R c ps 3.5) under condi- 
tions in which the Oberb e ck-Bo ussinesq (OB) approximation was quite well satisfied 
ijAssenheimer fc Steinberg (^996)). In that case hexagons with up-flow and those with 
down- flow in the center are equivalent and both were observed to coexist in adjacent 
domains. A subse quent numerical s tabilit y analysis confirmed the existence of such sta- 
ble OB-hexagons (|Clever fc Bussel l|l996h : iBusse et all i|l999l) ). They are characterized 
by a pronounced spatial concentration of the flow in the off-center regions. To resolve 
the strong spatial variations of the hexagons in the plane modes with relatively high 
wavenu mbers, which are negl ected in the familiar amplitude equations, need to be re- 
tained l|Busse fc Cleverl l)l999|) V 

In the other experiment, which used SFq near the thermodynamic critical point as the 
working fluid, it was found that the hexagons that arise in the immediate vicinity of the 
onset of convection can become unstable as the Rayleig h number is increased an d then 
restabilize again at higher Rayleigh numbers, e = 0(1) l|Rov fc Steinberd l)2002j) ). This 
restabilization was termed 'reentrance '. As the non-Boussincsq effects were increased the 
intermediate e-range over which rolls were the preferred planform shrank and eventually 
hexagons were found to dominate rolls from onset all the way to e = 0(1). Since the 
reentrant hexagons have been observed even for moderate NOB-cffccts the restabilization 
at larger e has been attributed by the authors to the high compressibility of SFq near its 
critical point ilRov fc Steinberg! 1120021) 1. 

In previous numerical stability calculations of rotating non-Boussinesq convection in 
water it was found that hexagons can be linearly stable over the whole range < e < 1 
ijYoung et all 1 2003)1. In the strongly nonlinear regime e = 0(1) a chaotic state ('whirling 
chaos') was obtained in which individual hexagonal cells oscillate or rotate, often inducing 
the nucleation of additional cells. Even in the presence of lateral walls, which typically 
induce the nucleation of rolls, this hexagon-based spatio-temporally chaotic state was 
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found to persist. In that investigation the mechanism that is responsible for the linear 
stability of the hexagons over a fairly wide range in e was, however, not understood. 

In the present paper we identify easily accessible parameter regimes of non-Boussinesq 
convection in which strongly nonlinear hexagons are linearly stable. We point out two 
mechanisms that contribute to their stability. First, the mechanism that stabilizes OB- 
hexagons at larger Raylcigh numbers also enhances the stability of the NOB-hexagons. 
Second, with increasing Rayleigh number the temperature difference across the fluid layer 
increases and with it the strength of the NOB-effects. A simple weakly nonlinear model 
shows that this effect alone can be sufficient to lead to a restabilization of the hexagons 
already at relatively low values of e, which experimentally leads to reentrant hexagons. 
Together, these two mechanism can lead to a restabilization for values o f e as low as 
e = 0. 2. Since our computations are based on water as a working fluid (cl . lYoung et all 
(2003)), which is essentially incompressible in the investigated regime, it is clear that 
high compressibility is not a necessary condition to obtain reentrant hexagons. 

The paper is organized as follows. In Sec|2]we briefly review the basic equations that 
we use, pointing out in which way our computations focus on weakly non-Boussinesq, 
but strongly nonlinear convection. In Sec|3|we present our results for the linear stability 
regimes of hexagonal and roll patterns. The two mechanisms leading to restabilization are 
discussed in SecQ] Direct simulations of the temporal evolution for boundary conditions 
mimicking circular containers are discussed in Sec|3] Conclusions follow in SecEl 



2 Basic Equations 

We consider a horizontal fluid layer of thickness d, density p, kinematic viscosity v, 
heat conductivity A, and specific heat c p . The layer is infinite in the horizontal extension 
and is limited in the vertical direction by two horizontal, rigid plates with high thermal 
conductivity. The system is heated from below (at temperature T\) and cooled from above 
(at temperature T2 < 7j ) . The governing equ ations expressing the balance of momentum, 
mass, and energy are ijChandrasekharl (§1,961) ) 

dt(pui) + dj(pujUi) = -d k p - pg5 i3 + dj (up (diUj + djUi)) (2.1) 
d t p + d j (pu j ) = (2.2) 

d t T + UjdjT = dj(XdjT). (2.3) 

pc p 

Here u = (7/1,1*2,1/3) is the fluid velocity, T the temperature, p the pressure, g the 
acceleration of gravity, and Sij is the Kronecker delta. Summation over repeated indices 
is implied. The origin of a Cartesian coordinate system with the z-axis perpendicular to 
the horizontal plates is fixed in the middle of the layer. As usually, viscous heating and 
volume viscosity effects can safely be neglected. 

Realistic rigid boundary conditions are taken at the two boundaries for the velocity, 

u = 0atz = ±- (2.4) 

and fixed values for the temperature 

T = T 1 =T + — atz = --, (2.5) 

T = T 2 = T - — at z = +-. (2.6) 
Here T = (T x + T 2 )/2 denotes the mean temperature and AT = (T x - T 2 ) > is the 
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temperature difference across the layer. We assume an experimental procedure in which 
To is kept constant while the main control parameter AT is varied. 

We focus in this work on weakly non-Boussincsq convection, so we keep the temperature 
dependence of the various fluid properties to leading o rder, and ex pand them about the 
mean temperature To in line with Busse's convention l|Bussei l)l967(l ') 

= l-7o—^ — (l+7i— — ) + •■■ (2-7) 

Po J s J s 

v(T) 



vo 
MT) 

Ao 
Cp(T) 

CpO 



1 + 7*-^ + ... (2.8) 
1+73^^ + ... (2.9) 
1+74^777^ + - (2-10) 



where po, fo, Ao, and c p o denote the values of the respective quantities at the mean 
temperature T). The dots denote higher-order terms to be neglected. Introducing the 
thermal diffusivity ko = \o/poc p , the scaling temperature T s = UoKo/(aogd 3 ) is used to 
define the nondimensionalized slopes of the density, viscosity, heat conductivity, and heat 
capacity at To in terms of the ji For instance, the usual heat expansion coefficient at 
T = T) is given by ao = 70 /Ts. Beyond the Boussinesq approximation also the curvature 
of p(T) at To, which is proportional to 7071/T, 2 , comes into play. 

To make the governing equations and boundary conditions dimensionless, the following 
scales are selected: for the length d 7 for the time cP/kq, for the pressure po^o^o/d 2 , and 
for the temperature T s . This gives rise to two dimensionless quantities: the Prandtl 
number Pr = vq/kq, an d the Rayleigh number R = AT/T S = a ATgd 3 / (voKo). Thus, 
we use in the following the dimensionless temperatures T = T/T s and T = T /T s , heat 
conductivity A = A/Ao, density p — p/po, kinematic viscosity v — v/vo, and specific 
heat c p = Cp/cpo- Finally, we write the equations in terms of the dimensionless momenta 
Vi = pu{dj poK-o instead of the velocities. In the following we omit the hats again for 
simplicity. 

Since the fluid ve locities a re sma ll compared to the sound velocity we make the anelas- 
tic approximation neglect the time derivative in the continuity 



equation (|2.2(l . This simplifies the computation considerably since it reduces the number 
of evolution equations. Furthermore, Vi becomes a solenoidal field, which can be repre- 
sented in the standard poloidal-toroidal decomposition by two velocity potentials l|Bussd 
(1989)) automatically enforcing the mass conservation. 

The (dimensionless) conduction solution (v = 0) of (Ol with 12. 512.612. 712.912. ll)f) is 
given by 

T cond = To + r( k -z- ^-(z 2 - 1) + 0( 73 2 )) • (2.11) 

We rewrite the temperature T in terms of the deviation & from the conductive profile 
neglecting its ©(7! )-contribution, 

9 = T - T cond = T - To - R (-z - y (z 2 - i)) . (2.12) 
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We then obtain as the final dimensionless equations 

(dm + v 3 d 3 {^\\ = -d iP + s i3 (1 + 71 (-2^ + 1) ) h 



P P 



(2.13) 



djVj = 0, (2.14) 

d t O + —d 3 @ = —d^XdjQ) - 7 3 <9 2 6 - R— (1 + j 3 z). (2.15) 
9 P c p P 



The dimensionless boundary conditions are 



v(x, y, z, t) — 9(x, y, z, t) = at z = ±— . (2-16) 
The nondimcnsionalizcd fluid parameters H2.7[) - I|2.10[l read now: 

p(9) = l- 7 o(-*+§), (2.17) 

K©) = l + 72(-^+|), (2.18) 

A(e)=l + 73 (-2+|), (2.19) 

c p (6) = l+74(-^+|). (2.20) 

where we have introduced the quantities 7; = 7^ = ^fiAT/T Sl which can be expressed 
as follows: 

Here 7 2 c = jiAT c /T s denotes th e non-Boussi nesq coefficients evaluated at the onset of 
convection, T = T c , as used in iBussd lll967|) . In our case the 7? = ^fiAT c /T s can be 
well approximated using the Boussinesq value for the critical temperature difference, 
AT C = 1708 T s . In principle, the non-Boussinesq effects induce small corrections through 
a shift in the threshold, which are easily numerically computed, if necessary. Note that in 
contrast to the 74 the 7, (AT) are linear in the main control parameter AT and therefore 
have to be adjusted in the nonlinear regime. 

We consider the non-Boussinesq effects to be weak and keep in all material properties 
only the leading-order temperature dependence beyond the Boussinesq approximation. 
Therefore the 71 -term appears explicitly in H2.13|l . while in all other terms it would 
constitute only a quadratic correction just like the terms omitted in H2.7(l - H2.10|l . Corre- 
spondingly, we expand the denominators in I|2.13I2.15|I that contain material properties 
to leading order in the 7*. In analogy to IBussd l|1967|) . we further omit non-Boussinesq 
terms that contain cubic nonlinearities in the amplitudes Vi or 0, as they arise from 
the expansion of the advection terms Vjdj(vi/ p) and (vj/p)dj& when the temperature- 
dependence of the density is taken into account. Since we will be considering Rayleigh 
numbers up to twice the critical value, which implies enhanced non-Boussinesq effects, 
these approximations may lead to quantitative differences compared to the fully non- 
Boussinesq system, even though the temperature-dependence of the material properties 
themselves may quite well be described by a linear (or quadratic in the case of the den- 
sity) approximation. Unfortunately we are not aware of rigorous solutions of Eqs (2) that 
would allow quantitative tests of our approximations scheme for finite ji. 
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In the weakly nonlinear regime hexagon and roll patterns are described by the ampli- 
tudes At of the three dominant Fourier modes associated with the wavevectors qi = 
g(l,0), q2 = (?/2(— 1, V3), Q3 = <z/2(— 1, — v3)- Roll solutions correspond to Ai ^ 
0, ^2 = A3 = and hexagons to A\ = A2 = A3 7^ 0, while for mixed hexagon so- 
lutions one has -A] ^Aaa^ 0. The amplit udes satisfy the well-known coupled amplitude 



with the equations for ^2,3 obtained from i|2.22|l by cyclic permutation. The calculation 
of the coefficients 5, gi involves vertical averages over certain products of the critical 
eigenvector components, which are obtained from a linearization of i|2.13l2.15|) . 

The quadratic coefficient S arises f r om th e NOB-effects. To leading order in the 7? 
it has been calculated first by IBussel l)l967(l and is proportional to Busse's parameter 
Q, which is conventionally used as a measure for the strength of the NOB-effects. It is 
defined as 



w here the Vj are certa i n line ar functions of Pr 1 , which can be found in the review article 



bv lBodenschatz et all l|2000r) t. Busse's parameter Q characterizes the breaking of the up- 
down symmetry, which renders at most one of the two possible types of hexagons stable. 
Gases have a positive value of Q and exhibit hexagons with down flow in the center (g- 
hexagons), whereas liquids have negative Q and show hexagons with up flow (Z-hexagons). 
In the weakly nonlinear approach the cubic coefficients gi in (|2.22|) are evaluated with 
-ff = 0, which is consistent with the assumption of small 7 4 c , 7? = O(Aj). In the case o f 
finite 7? the calculation can be refined along the lines presented bv lPlaut & PeschHl999|) . 

The stability of weakly nonlinear roll and hexagon patterns is determined by a linear 
stability analysis of the various solutions of Eq. I|2.22[l . The resulting bifu rcation diagram 
is sketched schematically in FigQ]for the case 52 > ffi fcf. IBussel (fl9 6711 . The hexagons 
arise unstably from the conductive state in a transcritical bifurcation and become stable 
through a saddle-node bifurcation. When the heating, i.e. e, is increased they become 
unstable in a transcritical bifurcation at en involving a mixed-mode solution. The rolls 
in turn are unstable at threshold and become stable through a pitchfork bifurcation at 
en, when the control parameter is increased. 

We focus in this paper on the stability properties of the patterns in the strongly 
nonlinear regime. They are deter mined by a Ga l erkin e xpansion of all fields in Eq. 
(I113I2M4I2MRI) (see, for instance. iBusse fc Cleverl (Il979h : IBussel il<M& ). Their depen- 
dence on the vertical coordinate z is captured by expanding them in appropriate com- 
binations of trigonometric a nd Chandrasekhar fun c tions in z th at satisfy the top and 
bottom boundary conditions l|Chandrasekharl l|l96lj) : lBussel jl967|) 1. In most of the com- 
putations we used n z = 6 modes for each field. With respect to the horizontal coordinates 
in the lateral directions we use a Fourier expansions on a hexagonal lattice. The Fourier 
wave vectors q are constructed as linear combinations of the hexagonal basis vectors 
bi = q(l,Q) and b2 = q(l/2, V3/2) as q = mbi + nb2 with the integers m and n in 
the range |mbi + nb2| < n q q. The largest wavenumbcr is then n q q and the number of 
Fourier modes retained is given by 1 + 6j2j=i 3- Typically we used n q = 3. Thus, solv- 

f The expressions given on p. 742 of Bodcnsch atz et al\ l2000f) use the symbol V instead of 
Q and correct a small error in Busse's calculation of Vs. 




(2.22) 
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Figure 1. Sketch of the bifurcation diagram for hexagons and rolls in the weakly non-linear 
regime of NOB-convection. The solid lines correspond to stable states and dashed lines to un- 
stable ones. The unstable mixed mode destabilizes the hexagons and stabilizes the rolls. 




ing the PDEs l|2. 1312. 1412. 15|1 is reduced to solving a system of nonlinear ODEs in time 
for the Galerkin expansion coefficients. The standard linear analysis of the ODEs yields 
the critical Rayleigh number R c as well as the critical wavenumber q c . Both depend on 
the NOB-coefficients 7? which in turn depend on R c . Thus, in principle one obtains an 
implicit equation for the 7?. The shift in the critical Rayleigh number away from the 
classical value R c — 1708 due to the NOB-effects is, however, quite small (less than 
1 percent) and therefore the resulting change in the 7? is negligible. In this paper we 
therefore choose the 7? corresponding to R c = 1708. 

To investigate the nonlinear hexagon solutions, we start with the standard weakly non- 
linear analysis to determine the coefficients of the coupled amplitude equations (|2.22() . 
To obtain the fully nonlinear solutions we need to solve the ODEs for the coefficients of 
the Galerkin expansion, which become a system of nonlinear algebraic equations in our 
stationary case. This is achieved with a Newton solver for which the weakly nonlinear 
solutions serve as convenient starting solutions. The solutions are tested for amplitude sta- 
bility by monitoring the growthrates of linear perturbations of the expansion coefficients. 
Since the system is spatially periodic (characterized by Fouriermodes with wavevectors q) 
the possibility of side-band instabilities with respect to modes with wavevectors q±s has 
to be considered as well. This is achieved by introducing Floquet multipliers exp(is- (x, y)) 
in the Fourier ansatz for the linear perturbations of the Galerkin solutions. 

We also study the dynamics of complex patterns that arise from instabilities of the 
periodic states. For that purpose we h a ve extended our previous ly developed spectral 
code for the OB-equations l|Peschl l|l99fil) : lBodenschatz et all l(2000l) ) to include the NOB- 
effects in Q2. 1312. 1412. 15JI . It employs the same vertical modes as the Galerkin stability 
code but places the wave vectors of the Fourier modes on a rectangular rather than a 
hexagonal grid. To solve for the time dependence we have chosen a fully implicit scheme 
for the linear terms, whereas the nonlinear parts are treated explicitly (second-order 
Adams-Bashforth method). The time step is typically taken to be i„/500, where t v is 
the vertical diffusion time. We have tested that the stability regimes obtained from the 
Galerkin analysis are consistent with the direct numerical simulations. 
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3 Linear Stability of Hexagons 

Instead of extensive parameter studies, we present in this work specific, interesting 
scenarios that should be experimentally realizable. We focus our investigation on wa- 
ter, which has a moderate Prandtl number and for which reentrant hexagons should be 
readily accessible in convection cells with conventional layer thickness d in a range of 
temperatures close to room temperature. 

3.1 Amplitude Instabilities 

In our analysis, we first concentrate on spatially periodic solutions with the wavenumber 
fixed at the critical wavenumber and discuss their domains of existence and stability as a 
function of the control parameter e = (R— i2 c (7?))/i? c (7?). We have chosen three different 
cells with thickness d — 1.5, 1.8, and 2.1mm, respectively. The case d = 1.8 mm is of 
partic ular interest since it has been studied in previous experiments l)Panrpaloni et al\ 

Tab. ngi ves the non-Boussinesq coefficients and the value of the non-Boussinesq pa- 
rameter Q at the onset of convection for a representative range of mean temperatures 
To in a fluid layer of thickness d — 1.8 mm]. As indicated before, the 7$ are linear in the 
temperature difference AT f see 12.21)) and therefore they depend on e, 

7i = 7° + =7f (1 + e). (3.1) 

Note that with increasing mean temperature the critical temperature difference AT C 
(given in the second column of Tab. ^) decreases. Therefore the variation of the fluid 
properties across the layer at the critical temperature and with it the coefficients also 
decrease with increasing temperature. Since we keep the mean temperature Tq fixed when 
changing the Rayleigh number, the accessible range in e is limited by the requirement 
that the temperature at the top plate, T(z = d/2) = T) — AT/2, be above freezing. It 
turns out that the full range < e < 1 is then only accessible for average temperatures 
T) above a certain temperature Tq , which is 20°C for d = 1.8 mm. 

Using the Galerkin method, we extend the weakly nonlinear result sketched in Fig^ 
ijBussd l|l967j) ') to the strongly nonlinear regime. Fig[2]shows the resulting stability limits 
for hexagons and for rolls. As predicted by weakly nonlinear theory, the hexagons are 
linearly stable with respect to amplitude perturbations for very small e with respect am- 
plitude perturbations. For not too small values of the mean temperature To and layer 
thickness d they become unstable as the control parameter is increased. This instability 
corresponds to the transcritical bifurcation at e# in Fig^ For the convection cells in- 
vestigated here, the hexagon patterns then undergo a second steady bifurcation as the 
control parameter is increased further and become stable again. As the mean tempera- 
ture is decreased or the layer thickness is decreased the critical heating and with it the 
non-Boussinesq effects increase. This shifts the point of restabilization to lower e and 
the lower stability limit to higher e, decreasing the e-range over which the hexagons are 
unstable, until the two limits merge at a temperature T m . For To < T m the hexagons are 
amplitude-stable over the whole range of e considered (0 < e < 1). 

The restabilization of hexagons in the strongly nonlinear regime has been observed by 
Roy and Steinberg i n S Fp, near the thermod ynamical critical point, where it has been 
termed 'reentrance' (|Rov fc Steinberg! |2£)02)). There it was argued that since the non- 
Boussinesq effects in that system are not very large near onset the reentrance is due to 

f These values were obtained with a code kindly provided by G. Ahlers. 
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12.94 


5.68 


0.0036 


0.2122 


-0.2725 


0.0352 


-0.0013 


-1.837 


32 


10.74 


5.18 


0.0034 


0.1440 


-0.2126 


0.0273 


-0.0005 


-1.292 


36 


9.12 


4.76 


0.0032 


0.1023 


-0.1709 


0.0216 


0.0001 


-0.954 


40 


7.87 


4.38 


0.0030 


0.0755 


-0.1405 


0.0173 


0.0004 


-0.731 


50 


5.75 


3.62 


0.0026 


0.0400 


-0.0926 


0.0104 


0.0007 


-0.428 


60 


4.42 


3.05 


0.0023 


0.0245 


-0.0654 


0.0064 


0.0006 


-0.287 



Table 1. Values of the Prandtl number Pr, non-Boussinesq coefficients "ff, and Busse's pa- 
rameter Q for water at the onset of convection as a function of the mean temperature and the 
temperature difference. The liquid layer has a depth of d = 1.8 mm. 



the large compressibility of the fluid in this parameter regime. By assuming that the 
working fluid is incompressible, which is an excellent approximation for water, our com- 
putations show that high compressibility is not needed for reentrance. The significance of 
the compressibility for th e occurrence of the reentrance has also been called into ques tion 
recently bv lAhlerd l)2005(l (see his footnote 59) based on the work bv lOh et al\ l|2004|) . 

We have also computed the stability of rolls with respect to amplitude perturbations. 
The corresponding stability limits are indicated in Fig[21 by thin lines. Note that it is 
not meaningful to extend these stability limits to lower values of To than shown in FigEl 
since then the temperature at the top plate would be below the freezing temperature of 
water. Below the thin curves rolls are unstable, but they become linearly stable when e is 
increased beyond the lines. This stability limit corresponds to the pitch-fork bifurcation 
at €r in Fig^ As the non-Boussinesq effects become stronger the stabilization of rolls is 
shifted to larger e. In contrast to the hexagons, for the convection cells investigated here 
the rolls do not undergo a second bifurcation, which would destabilize them, and remain 
amplitude-stable up to the largest values of e considered. For strong non-Boussinesq 
effects one has therefore a very large range of parameters over which the competing rolls 
and hexagons are both linearly amplitude-stable. 



The amplitude-stability limits of the hexagons and rolls depend, of course, on their 
wavenumber. This is illustrated in Fig. [31 where we fix the mean temperature T and 
determine the stability limits of the hexagons as a function of their wavenumber q. In- 
terestingly, for To = 36 °C (Fig 3(a) \ the hexagons become more stable for small and for 
large wavenumbers and the instability region forms a bubble-like, closed curve, outside 
of which the hexagons are stable with respect to amplitude perturbations. These bubbles 
are hyper-surfaces in q — e — Tb-space. With decreasing non-Boussinesq effects (increasing 
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Figure 2. Stability regions for water with respect to amplitude perturbations for three fluid 
depths: d = 2.1 mm (dotted lines), d = 1.8mm (full lines), d = 1.5mm (dot-dashed line). Thick 
curves: stability boundaries for hexagons. Thin curves: stability boundaries for rolls. For a given 
depth, rolls are stable above the thin line, and hexagons unstable in the inner region of the thick 
line. Stability limits obtained for the critical wavenumber q c . 



Tq) the bubble grows and extends to larger values of e. Eventually, the upper part of 



the bubble is shifted to e-values beyond the range considered in this paper (Fig 3(b) for 
T = 40°C). 

3.2 Side-Band Instabilities 

Using the Galerkin method, we have studied the stability of the hexagons with respect to 
long- and short-wave perturbations as shown in Fig|3] We find that over the whole range 
< e < 1 the only relevant side-band perturbations are long-wave and steady, as is the 
case in the weakly nonlinear regime. The long-wave perturbations involve longitudinal 
and transverse phase modes, which c a n be described by two coupl e d evo lution equations 



and transverse phase modes, which can be described by two coupled evolution equations 
l Lauzeral efali l|l993(k iHovlel (|l995h : lEchebarria fc Perez-Garcial |l998D. The stability 



limits obtained from the Galerkin analysis are shown in Fig. 01 I n this parameter regime 
the stability region consists of two disconnected domains, reflecting the reentrant nature 
of the hexagons. The stability domain near onset is very small and closes up as the 
amplitude stability limit is r eached. This b e havio r corresponds to that obtained from 
the weakly nonlinear theory l)Lauzeral et all l)l993h ). In the reentrant regime the stable 
domain opens up again in an analogous fashion when the amplitude-stability limit is 
surpassed. Note that the stability boundaries lean toward lower wavenumbers. Thus, 
stable reentrant hexagonal patterns are expected to have wavenumbers below q c . 

As the mean temperature is lowered the bubble of the amplitude instability shrinks 
and eventually the bubble disappears (Fig0J. The side-band stability limit reaches then 
without interruption from the strongly nonlinear regime all the way down to threshold 
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(b) 

Figure 3. Stability regions for hexagons in water with respect to amplitude (dashed line) and 
side-band perturbations (dot-dashed line). The depth of the fluid layer is d = 1.8 mm. Figure 
(a) To = 36 °C, (b) To = 40 °C. Hexagons are stable with respect to amplitude perturbations 
outside the dashed-line region, and stable with respect to side-band perturbation inside the 
dot-dashed line region. 
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Figure 4. Stability regions for hexagons with respect to side-band pertubations, in water. 
The fluid depth is d = 1.8 mm. Solid line: T = 28°C. Dashed line: T = 32°C. 

(more precisely to the saddle- node bifurcation of the hexagons). As the non-Boussinesq 
effects become yet stronger the range of stable wavenumbers widens. 

4 Origin of Reentrant Hexagons 

At first the appearance of stable reentrant non-Boussinesq hexagons seems quite sur- 
prising, in particular, considering the relatively small values of e for which the resta- 
bilization of the hexagons can occur. We have identified two major factors that con- 
tribute to their appearance. One factor is the fact that even in the Boussinesq case 
hexagons can be stable for sufficiently large Rayleigh numbers (above e ps 1). They have 
been observed in convection experiments using SFq close to its thermodynamical critical 
point as a working fluid where th ey nucleated in the cores of target and spiral patterns 
(jAssenheimer fc Steinberel (^996)). A subsequent numerical stability analysis confirmed 
the existen ce of stable OB-hexag o ns and attributed t heir appearance to the formation 
of plumes l|Clever & Bussel l|1996h : LBusse et all (|lj)99)). The second factor contributing 
to the reentrance of the non-Boussinsesq hexagons is the increase of the non-Boussinesq 
effects with the Rayleigh number (cf. (|3.1I) ). 

FigEl provides a quantitative assessment of the importance of the two mechanisms 
contributing to the reentrance of non-Boussinesq hexagons. Above the dashed-dotted line 
hexagons with wavenumber q c become amplitude-stable in the Boussinesq case 7,; = 0. 
Th is transition line corres ponds to the stability limit of Boussinesq hexagons found earlier 
hv lClever fc BuTsl lTl99^ (see their Fig.4). The very slight dependence of the stability 
limit on the mean temperature is due to the variation of the Prandtl number with Tq, 
which decreases from Pr = 5.4 at T a = 30 °C to P = 3.0 at T = 60 °C. With the 
non-Boussinesq effects included, the stabilization of the hexagons occurs at lower values 
of the control parameter. The dashed line in FigEl shows the resulting stability limits 
when the e-dependence of the 7* is neglected, 7$ = yf, while the solid line denotes the 
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Figure 5. Stability regions for hexagons with respect to amplitude perturbations in water for 
d = 1.8 mm. Solid line: stability limit when the dependence of 7 on the temperature is taken 
into account, 7; = 7 t c (l + e). Dashed line: stability limit when the 7* are fixed to their critical 
value, 7s = 7?. OB-hexagons are stable above the dashed-dotted line, 7^ = 0. Down-hexagons 
are stable above the dotted line for 7* = 7^(1 +e). 



stability limit with the dependence retained, ji = 7f(l + e). Thus, even when the non- 
Boussinesq effects are kept constant the stability limit connected with the stabilization of 
Boussinesq hexagons merges with the usual low-e stability limit when the non-Boussinesq 
effects become strong enough, i.e. at low temperatures. 

Since in the Boussinesq-case up- and down-hexagons are equivalent and become stable 
simultaneously it is to be expected that for weak non-Boussinesq effects both types of 
hexagons can become stable for large Rayleigh numbers, with the stabilization occurring, 
however, at different values of the Rayleigh number. For the up-hexagons, which in water 
are stable near onset, this stabilization corresponds to a reentrance. The stabilization 
of the down-hexagons is indicated in Fig|5] by a dotted line. As the non-Boussinesq 
effects become stronger the down-hexagons require ever higher Rayleigh numbers for 
stabilization. 

Even though the e-dependence of the non-Boussinesq effects is not the central driving 
force for the reentrance, it is instructive to discuss its effect on the hexagons briefly within 
a weakly nonlinear framework. We consider as a minimal model a slightly generalized 
version of Ea. (|2.22|) in which the quadratic coupling coefficient grows linearly with e, 

d t At = e 2 (ni • V?A 1 + eAt -(5 + fie)A 2 A 3 ~ 9l \A X \ 2 A 1 (4.1) 
- 52 (|A 2 | 2 + |A 3 | 2 )Ai 

Compared to (|2.22[) we also include a spatial gradient term involving the normal deriva- 
tive rii • V with x\.i — q,-/|qj|, i = 1, 2, 3. It allows long-wave modulations of the amplitude 
that capture side-band instabilities. Note that such a linear correction of the quadratic 
resonance term has been considered pr eviously in the context of hexagonal patterns in 
ferrofluids exposed to a magnetic field IjFriedrichs & Engell lj200lJ )). 
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Figure 6. Sketch of the bifurcation diagram obtained from 14, H (cf. Fig|7| below for k = 0). 



To get a complete picture of the impact of the e-dependence of the non-Boussinesq 
effects within the framework of the simple model (|4.1|l we introduce a rescaled amplitude 
A = g\A/{8 + fxe), a rescaled time t = (5 + iie) 2 t/gi, and a rescaled space variable 
f = (5 + iie)v / t;J~g~i. Eg. 1)4. II) can then be written as 



8 t A x = (n a • V) 2 A + eVU - A 2 A 3 - \A 1 \ 2 A 1 - ^ (\A\ 2 2 + \A\ 2 3 ) A 1 
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Here the rescaled control parameter is given by 



(6 + /ie) 2 



(4.2) 



(4.3) 



which implies that each value of e corresponds to two values of e, ei.2 = {gi — 25(ie ± 
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L (5^e)/2^ 2 e. 



For the case g\ > and /iS > 0, which is of interest here, a monotonic increase in e 
is mapped into a non-monotonic change of e with the range < e < e m = d / fi mapped 
onto < e < gi/45fj, and the range e m < e < oo mapped in a reverse fashion onto 
the same interval, gi/45/j, > e > 0. Thus, with increasing e the standard bifurcation 
diagram sketched in Figgis traversed towards the right up to e = e m . Note that the 
parameter e in Fig plays the role of e in ^4.2(1 . As e is increased further e decreases 
implying that the path through the bifurcation diagram in Figgis reversed. Thus, if e m 
is not too small the same mixed-mode that is created in the bifurcation stabilizing the 
rolls and that destabilizes the hexagons at a larger value of e in a transcritical bifurcation 
restabilizes the hexagons in a second transcritical bifurcation and eventually disappears at 
the bifurcation that destabilizes the rolls again. This scenario is shown in the qualitative 
bifurcation diagram depicted in FigO For larger values of fi, 



^ 1 {91-92? 

» > ^ = TS2(2g 1+ g 2 ) 



(4.4) 



e does not reach the stability limit of the hexagons and the hexagons remain linearly 
stable for all values of e. Specifically, for [i = /i m the two transcritical bifurcations at e# 
and e' H coincide and the upper and lower stability limits in FigEl merge, eliminating the 
amplitude- unstable regime of the hexagons. 
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Figure 7. Amplitude and side-band stability limits based on the extended 
Ginzburg-Landau-equation model 14. H with £ = 2.83, S = 0.048, [i = 0.28, g\ = 1, and 
g2 — 1.45. Shown are the neutral curve (dashed-dotted), the long- wave side-band stability limit 
(solid and dotted), and the amplitude stability limit (dashed). Hexagons are stable inside the 
solid lines. 



Fig. Ogives the side-band stability limits obtained within the minimal model 1)4.1(1 for 
a typical case. The solid lines denote the transverse long-wave phase mode, while the 
dotted line marks the longitudinal phase mode. The bifurcation diagram shown in FigEl 
corresponds to traversing this phase diagram at q = q c . Note that the destabilization 
of rolls at e^j is beyond the range of e shown in Fig. [7| In principle, either of the phase 
modes can lead to an instability. Within eq. I|4.ip . however, the longitudinal phase mode is 
relevant only for an extremely small range of e near threshold. For small Prandtl numbers 
eq. (14.1(1 would have to be extended to include a mean flow. Then the longitudinal phase 
mode can dominate the transverse mode over significa nt portions of the stability limits 
l|Young fc Rieckd (|2002|) ; ISemwogerere fc Schatzl (|2004|) ). 

Thus, the minimal model ((4.1(1 captures qualitatively the restabilization of the hexagons, 
the merging of the two stability limits as the non-Boussinesq effects are increased, and 
the main features of the side-band instabilities obtained in the full stability computations 
described in SecOl We have not found, however, the destabilization of rolls by a steady 
mode as it is suggested by the minimal model (cf. FigEJ). Instead, we find at quite large 
e an oscillatory instability. 

The full numerical stability analysis presented in SecOHdisplayed a tendency of the side- 
band-stable regions to shift to lower wavenumbers as e is increased. In principle, this could 
be mo deled phenomenologically by retaining nonlinear gradient te r ms in the minim al 
model ((Bragard & Velarde! l|l998l) : lEchebarria ^ fc Perez-rxarcfal|l99l : lNuz et al\$l99$) ). 
We will not pursue this here. Instead we point out that in experiments and in numerical 
stability analyses the Boussinesq hexagons are found to be st able with respect to side- 
band i nstabilities only for wav enumbers noticeably below q c (jAssenheimer fc Steinberg! 
(1996); IClever fc Bussd ()l996|) ). Considering the significance of the mechanism under- 
lying the stability of the Boussinesq hexagons for the reentrance of the non-Boussincsq 
hexagons, it is to be expected that the reduced wavenumber is a characteristic feature of 
this mechanism of reentrance. 
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Figure 8. Succession of snapshots for To = 24 °C in a circular cell of water of thickness 
d — 1.8 mm. The diameter of the cell is L = 8 ■ 2n/q c and the snapshots correspond to an 
integration time of 400 t v . 

5 Numerical Simulations 

To make closer contact with the results that would be expected in experimental inves- 
tigations we perform also direct numerical simulations of I|2.13I2.15[) . While our Galerkin 
approach uses realistic boundary conditions at the top and bottom plate it employs pe- 
riodic boundary in the lateral directions. Thus, these computations are able to predict 
instabilities that arise in the interior of the system, but they do not capture phenomena 
associated with the lateral walls. In most of the recent experiments circular containers 
have been used and typically it has been found that the walls prefer rolls, which are 
predominantly either perpendicular or parallel to the wall depending on details of the 
side-wall conditions. Thus, even very close to onset, where in the non-Boussinesq case 
rolls are unstable to hexagons, in these experiments a narrow ring of roll- like structures 
arises. 

A consequence of the wall-induced preference of rolls over hexagons is the experi- 
mental observation that the transition from hexagons to rolls does not occur at the 
Rayleigh number a t which the hexagons beco me linearly unstable, but already at lower 
Rayleigh numbers ijBodenschatz et cd\ 1^)91)). The transition is therefore more appro- 
priately interpreted as arising from a competition between rolls and hexagons, which are 
simultaneously linearly stable in that regime. The lowest-order Ginzburg-Landau equa- 
tion l|4.1|l is variational and therefore, within this framework, the transition is expected 
to occur when the energy of the rolls becomes lower than that of the hexagons. The full 
Navier-Stokes equations are not variational and therefore further away from threshold it 
is more appropriate to discuss the transition in terms of the invasion of one state into 
the other with the velocity of the front separating the two states going through zero at 
the transition point. 
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Figure 9. Results of simulations for a circular cell of thickness d = 1.8 mm and diameter 
L = 8 • 2 7r/q c . The simulations have been carried out for T = 20 °C, 24 °C, 28 °C. Circles 
correspond to rolls and triangles to hexagons. The stability limits for amplitude instabilities for 
hexagons (full line) and for rolls (dashed line) are reproduced from FigEl 



In order to mimic the experimentally employed circular containers we apply a strong 
radial subcritical ramp in the Rayleigh number that suppresses any convection outside 
a certain radi us. Such a bo u ndary condition is found to prefer rolls p erpendicular to 
this 'wall' (cf. iDecker et all ljl994h . Fig.l in iBodenschatz et all (|200nh ). Snapshots of 
simulations in such a circular cell with diameter L = 8 ■ 2n/q c and layer thickness 
d = 1.8 mm are shown in Fig. [H] In all cases random initial conditions were used. Near 
threshold (e = 0.1) the cell is completely filled with regular hexagons, except for a 
narrow ring of roll- like structures that are driven by the boundaries. When the heating is 
increased (e = 0.2) rolls in the interior become linearly amplitude-stable (cf. FigEJ) and 
the rolls that are driven by the boundary invade most of the c ell. So far the scenario is 
very similar to the experimental observations shown in Fig. 4 of IPampaloni et al .] lll992l) 
for To = 28°C (see also Fig.[5]below). While the roll pattern persists near the boundaries 
for yet larger values of e, a domain of ordered hexagons appears in the center of the 
cell for e = 0(1). The size of the inner domain of hexagons grows with increasing e and 
for e = 1.4 the hexagons fill essentially the whole convection cell. The growing of the 
hexagon domain with increasing e can be understood to arise from a balance between the 
increasing tendency of the hexagon domains to invade domains of rolls on the one hand 
and the predominance of rolls near the boundaries on the o ther hand. It is worth noting 
that according to our computations IPampaloni et all (|l992j) should have obtained stable 
reentrant hexagons if they had gone to larger values of the Rayleigh number. 



A comparison of the snapshot at e = 0.1 with that at e = 1.4 shows that the reentrant 
hexagons arising from random initial conditions have a smaller wavenumber than the 
hexagons that appear at threshold. This trend is consistent with the results of the side- 
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band stability calculations (Figs.ElE}, which show that with increasing e the wavenumbcr 
range of stable hexagons moves toward lower wavenumbers. 

In Fig.|^|the results of a set of simulations in a circular cell with the same dimensions 
as in Fig. are summarized for three mean temperatures, To = 20 °C, Tq = 24 °C, and 
To = 28 °C. Again random initial conditions are used. The triangles denote parameter 
values for which the final state consists mostly of regular hexagons, while the circles 
indicate a final roll state. When both patterns coexist over the course of the simulation 
(tmax = 200i„) both symbols are plotted. The simulations show that for To = 20 °C the 
NOB-cffccts arc so strong that hexagons dominate rolls over the full range of e studied. 
As seen before in Fig|51 for a mean temperature of Tq = 24 °C an intermediate range of e 
arises in which the final state consists of rolls. For e = 1.0 rolls and hexagons coexist and 
for yet larger e reentrant hexagons appear. If the mean temperature is increased further 
to T) = 28 °C the same sequence: hexagons — * rolls+hexagons — * rolls — > rolls+hexagons 
— > hexagons is obtained. The only difference to the case To = 24°C is that the reentrance 
is shifted to larger values of e. 

To make contact with the linear stability analysis for periodic boundary conditions 
Fig. shows also the stability limits for hexagons and rolls with solid and dashed lines, 
respectively. For the comparison with the results of the simulations in the circular con- 
tainer it is important to keep in mind that the point at which the hexagons and rolls have 
equal energy, or more generally, at which the velocity of fronts separating the two states 
changes sign is not given by the linear stability limit of either the hexagons or the rolls, 
but rather lies somewhere between the two stability limits. It is therefore reasonable to 
expect that as the upper linear stability limit of hexagons is shifted to lower values of the 
Rayleigh number with decreasing mean temperatures To, the Rayleigh number at which 
hexagons start to invade rolls also decreases and correspondingly the transition from rolls 
to hexagons in the presence of boundaries is also shifted to lower Rayleigh numbers. 



6 Conclusions 

In this paper we have studied non-Boussinesq convection in water for realistic param- 
eters and boundary conditions. We have complemented numerical stability analyses for 
periodic boundary conditions with direct numerical simulations that mimic set-ups used 
in usual laboratory experiments. Our main result is the finding of reentrant hexagons, 
i.e. we fin d that the hex agon patterns, which typically become unstable not too far from 
threshold can regain stability further above threshold via a specific resta- 

bilization transition. For strong, but realistic NOB-effects hexagons that become unstable 
at e = 0.15 can become stable again already at e = 0.2. For yet stronger NOB-effects 
hexagons are amplitude-stable over the whole range of e investigated (e < 1.5). This 
stabilization over the whole range of e is not due to a shifting of the initial transition 
from hexagons to rolls to ever increasing e, but rather to the collision of this transition 
with the restabilization transition of the hexagons. 

Reentrant non-Boussinesq hexagons have been observed in convec tion experiments 
using SFq near its thermodynamic critical point as a working fluid l|Rov fc Steinberg! 
(2002)). Their restabilization has been attributed to the strong compressibility in this 
regime. Our computations, being based on water as the working fluid, demonstrate that 
compressibility is not necessary for this phenomenon. We show that in water the reen- 
trance is instead connected with the f act that even in the Bouss i nesq-c a se hexagons can 
becom e stable for sufficiently large e ((Assenheimer fc Steinberg! l|l99fi(k IClever fc Bussel 
l|l99fil) ) and with the increase of the non-Boussinesq effects with e. A simple amplitude- 
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equation model capturing the latter aspect provides qualitative insight into the stability 
of the hexagons, including their side-band instabilities. 

Reflecting the fact that Boussinesq he xagons are stab l e with respect to side-band 
perturbations only for low wavenumbers l|Clever fc Bussel <|199^ V we find that in wa- 
ter the wavenumber of the reentrant hexagons is noticeably below q c . In contrast, the 
strongly nonlinea r hexa gons found experimentally in SFq do not exhibit this trend 
dRov fc Steinberg! ll2002ln . Since we find the same tende ncy towards lower wavenumbers 
also in gases away from the thermodynamic critical point ijMadruea et al\ l|in preparation^ 
one might speculate that the stabilization of the experimentally observed hexagons is 
due to a mechanism that differs from the mechanism of reentrance discussed here. In 
the experimental S'Fg-system even small changes in the mean density can modify the 
non-Boussinesq effe cts sign i ficant ly, which would strongly affect the transitions between 
rolls and hexagons <|Ahlersl tOOti ). Compared to SFq near the critical point, water has 
the great advantage that the non-Boussinesq effects are much easier to control. 

The connection between the reentrant hexagons and the Boussinesq hexagons suggests 
that even in the non-Boussinesq case down-hexagons may become linearly amplitude- 
stable, albeit for yet higher Rayleigh numbers. We show that this is indeed the case. We 
have not investigated their side-band instabilities and it is not clear whether they can 
coexist with up-hexagons or whether domains of up-hexagons always invade domains of 
down- hexagons. 

To address the stability and dynamics of the hexagon patterns in large-aspect ratio 
systems with non-periodic boundary conditions we have also performed direct numerical 
simulations of the Navier-Stokes equations. By a suitable strong variation of the local 
Raylei gh number we have i mplemented a qualitatively convincing model of a circular con- 
tainer ( Decker et al For inter mediate NOB-effects we c onfirm the experimentally 
observed scenario in the low-e regime l|Pampaloni et al\ l|l992() '). Our computations sug- 
gest that even in this experimental set-up reentrant hexagons should have been accessible 
experimentally for somewhat stronger heating (e « 1.5). 

Since in experiments the side walls always induce roll convection the transition be- 
tween rolls and hexagons occurs when the velocity of a front separating the two states 
changes sign and not at the stability limit of the hexagons. It is therefore not possible 
to study the amplitude-stability limit of the hexagons in such cells. We expect, how- 
ever, that by forcin g hexagonal patterns near the walls either by an appropriate space- 
depen dent heating ( Scmwogcrcrc fc Schatd ll20r)2h 1 or by a suitably corrugated bottom 
plate l|Bodenschat d ijunpublishedr) ') hexagons can be stabilized near the walls and the 
bulk stability of hexagons with respect to rolls can be investigated. 

The merging of the lower stability limit of hexagons with the restabilization line pro- 
vides also an explanation for the large contiguo us stability r a nge t hat was found in 
rotating non-Boussinesq convection using water ijYoung et al\ (|Jo03)). An interesting 
question is how the restabilization interacts with the Ho pf bifurcatio n of the h e xagon s 
to oscillating hexagons, w h ich is induced by rotation (ISwiftl <|1984l) : ISowardl l|1985h : 
lEchebarria fc Rieckel jiflfflll) : iMadruga fc Perez-Garcfal l|2004la 

We thank G. Ahlers for providing us with the code to determine the non-Boussinesq 
coefficient and for stimulating discussions. We gratefully acknowledge support by the 
office of Basic Energy Sciences at the Department of Energy (DE-FG02-92ER14303) and 
an NSF-SCREMS equipment grant (DMS-0322807). 
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